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Introduction 

It is well known that the presence of a real atomic level which is nearly resonant 
with the pump field can greatly enhance the Raman emission cross-section. In order to 
accurately calculate the Raman gain in systems where resonance enhancement plays a 
dominant role, expressions for the pump and signal susceptibilities must be derived. 
These expressions should be valid for arbitrary field strengths in order to allow for pump 
and signal saturation. In addition, the theory should allow for arbitrary longitudinal 
and transverse relaxation rates. This latter point is extremely vital for three level 
atomic systems such as the alkali earth metals since they do not have population reser- 
voirs and can have widely varying spontaneous lifetimes on the three pertinent transi- 
tions. Moreover, the dephasing rates are strong functions of electronic states and are 
therefore also different for the three coupled pairs of levels. These considerations are 
not as important when molecular systems are concerned since the large reservoir of rota- 
tional states serve to produce essentially equal longitudinal recovery rates for the popu- 
lation of the three levels. 

The most general solution to date is that of Temkin and Panock^ 1,2 ) who have 
solved the semiclassical rotating wave limit of a three level system under the constraint 
of a single longitudinal and transverse relaxation rate for the entire system of levels. 

We have for the first time solved the three level system with three arbitrary longi- 
tudinal and transverse relaxation rates. There is no need for setting either pair of rates 
equal and the expressions are valid for arbitrarily strong fields. The next step will be to 
velocity integrate these expressions for the cases of copropagating and counter- 
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propagating pump and signal fields in order to study the cesium system pumped by the 
doubled alexandrite laser emissions. 



Density Matrix Solutions for the Susceptibilities 
of a Three-Level System with Arbitrary Relaxation 
Rates and Field Strengths 
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The density matrix formalism is fully equivalent to the Schrodinger wave function however 
it is more applicable to the statistical case, in which the wave function is not known exactly. In 
this case the elements of the density matrix are taken to be ensemble averages of possible 
configurations of the system. As a result every wave function can be expressed as a unique den- 
sity matrix, but not all density matrices can be expressed as wave functions, reflecting the fact 
that the wave function contains only quantum uncertainties where as the density matrix con- 
tains statistical as well as quantum uncertainties. 

For a given wave function 

l*> = £C n |n> 


the density matrix is written as 

Py j = c7c; 


where the bar denotes the ensemble average. The time dependence of the elements is given by 

Pij = j^[p, H liJ 

where H is the total Hamiltonian of the system used to describe the three and four level atoms 
pictured below 



Figure 1 

Both configurations approximate electronic states in an alkali atom. The lower states jl> 
and |2> are s-states with principle quantum numbers n and n+1. In the four level atom |3> 
and |4> are p-states with j = 1/2 and 3/2 respectively and principle quantum number n. In the 
three level atom the splitting of the upper p-!evel has been neglected to simplify the calculations 
and results. 

The Hamiltonian in the dipole approximation is given by: 

H = H 0 - 7 z E! 

where H 0 is the unperturbed atomic Hamiltonian and -/ r® is the atom field Hamiltonian in the 
dipole approximation. We will assume that the dipole matrix element /j 12 is zero because both 
are s-states. 
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The time evolution of the density matrix elements for the three level system are therefore 
described by 

” ( j t ' — ■j~(Pi3 _ _ '7ii(/ , n _ Pn) 


dp22 

dt 


~ ^ h”^ 23 ~~ ^2)^23 '®( l ) ~ ^22(^22 ~ P22) 

— T" [(Pl3 ~ P3l)A*12 + (P32 _ /*23)/ i 23 J 'E'(t)'Y33(P33 _ P 33 ) 


Pl3 — k" £(P33 _ Pll)/*13 _ /*23Pl2 j ~ '(^31 ~ *713)^13 

P32 — j^(P33 ~ ^ 22)^23 ~ PnP\2 j 'E(t) + i(tU 32 - 1132)^32 

Pi 2 = Y [a‘13P32 - /‘ 32 P 13 ] ^(0 - »( w 21 “ i-7l2)Pl2 

where — (Ej - Ej)/h, E(t) — E p cos(cu p t 4- il> v ) + E 3 cos(a> 3 t + rp 3 ) and 7 ^ are phenomenological 
relaxation and dephasing rates, and p 3l , p 2 3 and p 12 are complex conjugates of /> 13 , p 32 and p l2 
respectively. The equations describing the four level system are given by: 

Pn — ^i3(P3i - P 13 ) + PwiPn ~ P 14 ) j 'E(t) - 7nPn 

P 22 ~ ^2o(P32 ~ P 23 ) + / i 2l(P42 _ P 24 ) j 'E(t) - 722 P 22 

P 33 = 'j~ ^13(Pl3 “ P 31 ) + A 1 23(P23 _ P 32 ) + p2i(P*2 ~ P 24 )j 'E(t) - 7 33 /? 33 

P44 = Y [^loCPH - P4l) + / i 2l(P24 ~ P 42 ) + / J 34(P34 ~ P 43 ) j 'E(t) - 7 44 P 44 

Pl2 = |a 1 13P32 ~ / i 23P33 "b M 14 P 42 - M 42 P 34 j 'E(t) - i(w 21 - ■ - 7l2)Pl2 

Pl3 — £/*13(P33 _ Pi l) + Ml4P43 ~ P23P12 ~ M43P14 j 'E(t) - i(w 3 i - i7is)pl3 

Pl4 = "j^ ^14^44 ~ Pll) + M13P34 _ M24P12 _ / 1 34Pl3 j 'E(t) - i(w 41 - i7l4)Pl4 

P 23 = ^M23(P33 _ P 22 ) + A*24P42 ~ P'U,P2l ~ l i 2iPi2 j 'E(t) - i(tU 32 ~ i^23)P23 

P24 ~ ^ j^/ i 24(P44 _ P 22 ) + / 1 23P34 ~ t l liP2 1 _ / i 34P23 j 'E(t) - i(utj 2 ~ 1724^24 
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Pii — ^ ^34(^44 £33) + / i 23/ , 14 ~ A i 23/ , 24 ~ PuPiX _ ^24^32 j - ^(t) ~ *( w 43 _ *734)^34 

It is obvious that these equations with their complex conjugates are much to complicated to 
solve. The remainder of this paper will be devoted to the three level system. 


Steady State Solutions 


To obtain the steady state solutions to the density matrix equation we first take out the 
explicit oscillations of the off diagonal elements 

P32 — Xe 


P 13 


Ae IWpt 


P12 = De 


'(<q> - ‘‘O 4 


and then neglect all terms which oscillate faster than u> p or u) s in the off diagonal element equa- 
tions and all DC terms in the diagonal element equations. This is the rotating wave approxima- 
tion which therefore disallows the possibility of Bloch-Seigert shifts in the resonances. If we 
then set all time derivatives equal to zero and define the following: 


fi,= 


2 h 




M13E, 


'V> D 


2 h 


L s = w s - w 32 + 1/^ = 6, + i / t x 
L p = w p - W31 - >/ r 2 = - i/r 2 

L ap == CJ 3 - Wp + W 2 1 + i/73 = ^ + i/r . 3 

the equations become 


L S X = &*A 3 2 - /? P * D 
L p A = -~/? p A 13 - 
LapD = ~/? p X + &*A 


and 


0 = 2 Im/? 3 X - 4 Im/? p * A + \ -J" + 

05 A 1 A 3 


(A 13 - A,° 3 ) + — — 



A32) 


0 = 2Im/?p A - 4Im/9 3 X + — 

3 



2 

T 2 


(A 32 - A32) + 


3 




We can now eliminate D from the first three equations and solve for X and A 

p : 


X = 


A 32 + |/?p| 2 R 


A ] 3 ^ A 32 



6 


^{a 13 - |&| 2 r 


A ] 3 + A 32 


Lp ^3 


where 


R = 


LgLp 


4 ft 2 - L p /?p 2 + L 3 L p L sp 

In terms of two real quantities, R, and R 2 , we have that: 


R = R, - iR 2 = 


(ip 2 + ^ 2 )(i s 2 + rf 2 )(A - iB) 
A 2 + B 2 


A = i 3 (i 2 + r 2 2 )(* s 2 + rf 2 - tf) - 5 p (i s 2 + rf 2 )(5 p 2 + rf 2 - tf) 

B = t;\6? + rf 2 )(* 2 + r 2 2 ) + + r 2 2 )(5 s 2 + rf 2 ) + fifaW + rf 2 ) 


This then yeilds 


and 


where 


Im/? 3 \ = 


Im/L*A 


IU 2 

-|/>T 

I L o I 2 


rf 1 + '*'* 


r.> - 


i Ls I 2 

1 1 A 11 
I L P I 2 


a 3 


°1 


\Pp\~\P,\ 

I Lg | 2 | L p | 


A I 1 2 1 A' 1 2 

I L. I 2 I L p | 2 


-a 2 A 


13 


a 2 A 


32 


«1 = 2 r 2 ‘^R, - ( 5 p - r 2 ~ 2 ) R 2 
q 2 = (Vi ~ <5 3 r o“ 1 ) R i + (n V + i p i 3 )R 2 
a 3 = 27! ^jRi - (5j 2 - r, 2 )R 2 

Substituting these into the last two density matrix equations then gives 


2 

1 

A 0 

1 

1 

1 

It, 

T 

A 3 

Ai 3 + 

3 

[t, 

" T 3 


2 1 ^3 1 2 

• _ ,0 

r* 1 + 2 | /? p | 2 

1 4 | * 

- 


A 0 
^32 


-Cto Q3 

“ + 


i L P I 2 


14 


+ 


4 1 h 1 2 
I L p f 


To 


l^l 2 


«1 q 2 


I L P I 2 


2 I L, | 


+ 


3 T, 


— f — 

3 It, 


1 

T 3 


*32 


+ 


l 13 
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H-L + JL A o + J_| A o 

3 T 3 T 2 32 3 I T 2 T 3 13 


4 ! ft I 


rf 1 + J ft 1 2 J 


~ a 2 Q3 

ft 1 2 ' Iftl 2 




2 1 1 2 
I L P I 2 


pi L-l I a I 2f Q 1 , - a 2 


r 2~ - I ft I 


2 1 |Lp| 2 I ft I 2 




Then by making the following substitutions and definitions 


W 2 , 1 ) , 4 I 1 2 ft-i , a2|2 [ «i 

3 ft t7) + tfTlFT 2 ■ lft 1 [«J77 


if 2 2(«, 2 +rf 2 ) 


, _ 2 lftl' ft- 1 , 0 , d | 2 T ~ ft 2 , Ilf 1 _ 1 

2 ft 2 + rf 2 \ 1 Mpl Sp + r 2 “ 2 2(ft + rf 2 ) i 3 ( T x T 3 

P “Iftl f _j I o 1 2 [ a i . 

r 3 — T7~ 77 { r 2 - I ft I T7— 77 + 


2 I ft, I 2 


ft + r 2 




, I 2 ft l ‘^ a 2 LlfJ L 

ft 2 + r 2 2 ft 2 + rf 2 J 3 l T 2 T 3 


r <ftft + 2 ft + 77^ft + 


-2\ c 2 1 -2' 

2 ) ft + r l 


The equations become 


— — — I — 

3 T x T 3 


+ 7 F - A 13 + 77 f 77 77T~ I ^32 — ftA 13 - r 2 A 32 


3 1 T, T 3 


j \ y ~ + y-|a 3 ° 2 + v(^-- 7^]^ 0 3 = -r 3 Ai 3 + r 4 A 32 

3 2 


3 [ T 2 T 3 


with the solutions 


Ai-j — 


r 4 \- [— + — )a, 0 3 + -(— - t~— ] Ag^] + r 2 [- 1— + — 

4 3 T t t 3 ) 13 3 ( t x t 3 j 32 J 2 3 ( t 3 t 2 

ftft - r 2 r 3 


3 ( ft T 3 


a 3 ° + — — 

3 “ 3 I T 2 


A, 1 


ft K4- + 4-|a 3 ° 2 + if^--^-]A 1 ° 3 +r 3 I(-±- + + 


3 ( ft T 2 


3 l T 2 T 3 


3 [ Tj T 3 


3 T 


ftft - r 2 r 3 


The polarization of the system can now be determined from the expectation value of the 
dipole operators by: 
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<p> = Tr(/tp) 

= 2Re(/i 31 p 13 ) + aRe(/i 32 p 23 ) 
which, after substituting for p 13 and p 32 becomes 


P = 2Re 


+ 2Re 


isr{ A - +l ^ | 2 R '[if + if]K"' 


i(w p t +i>) 


By comparing this equation with 


P + Re (x p E p e i( " pt+ ^) + Re ( Xs E s e iKt+W ) 


the complex susceptibilities are found to be given by: 

- I M13 I 2 ( . , . , 

Xp | Aj 3 I ft 3 I 


•2 D A 13 , A 32 

R V ~ 


+ I ^23 I 


N 


A 32 + l^rl 2 R* 


2r‘ II 


Lp L, j / 


In their complex forms 

X P = X P ' + ix P ' X 3 = Xs' + iXa' 


I Mis I 2 f [-1 l&| 2 «i 

’A,, - 


h(C+r 2 - 2 ) 


2 A _-2-\ I 13 2 (t 2 , _-2\ 


1 ^ . I ft 1 2 °2 
(S! + t 2 -*) 32 (S?+t{ 2 ) 


„ ,1 I /*32 I 2 

Xs ,/r 2 -2’ 


I ^ | r -l | I I °3 

|Ak| ‘ + WTTp) 


Xp = 


Xg 


| /x 13 1 2 r 

: 131 a 13 <5 p - 

h(<5p -P r 2 2 ) I 


I /? P I 2 «2 

k 13 , f o _Ov 

(*p + V) 


I ft 1 \ I ft 1 2 72 

(ft 2 + r 2~ 2 ) 32 (ft + T{ 2 ) 


I M32 I 
h(ft 4- i 




A32 ^3 


ft I H 


(ft + rf 2 ) (ft + rft) 


I ft I 2 7 2 

13 /c2 , _-2\ 


where 


7i = (ft 2 - ^ 2 )R! + 2r 2 ftR 2 
72 = (ftft + »fV)Ri + {5 s t 2 1 - i p r 1 “ 1 )R 2 
7 3 = (ft 2 - rf 2 )Ri - 2rf 1 ftR 2 

These susceptibilities are the first results which allow for completely arbitrary relaxation 
and dephasing rates as well as arbitrarily strong pump and signal fields. 
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